(12) United States Patent 

Heavlin 



iiiiiiiiiiiiiniiiiiiiiiiiiiiiiiiiiiiilHi 

US006366822B1 

(10) Patent No.: US 6,366,822 Bl 
(45) Date of Patent: Apr. 2, 2002 



STATISTICAL PROCESS WINDOW DESIGN 
METHODOLOGY 



(54) 

(75) Inventor: 
(73) 



William D. Heavlin, El Granada, CA 
(US) 



Assignee: Advanced Micro Devices, Inc., 
Sunnyvale, CA (US) 



( * ) Notice: Subject to any disclaimer, the term of this 
patent is extended or adjusted under 3^ 
U.S.C. 154(b) by 0 days. 

(21) AppL No.: 09/130,528 

(22) Filed: Aug. 4, 1998 

(51) Int. Cl7 G05B 13/02 

(52) U.S. CI 700/31; 700/121 

(58) Field of Search 700/121, 103, 

700/109, 31, 117, 29; 716/4 

(56) References Cited 

U.S. PATENT DOCUMENTS 

5,301,118 A ♦ 4/1994 Heck et al 700/109 

5,319,564 A * 6/1994 Smayling et al 703/13 

5,418,974 A * 5/1995 Craft et al 716/2 

5,526,293 A * 6/1996 Mozumder et al 716/19 

5,646,870 A * 7/1997 Krivokapic et al 716/4 

5,655,110 A * 8/1997 Krivokapic et al 716/19 

5,719,796 A * 2/1998 Chen 700/103 

5,787,269 A . * 7/1998 Hyodo 703/13 

OTHER PUBLICATIONS 

M. D. McKay and R. J. Beckman, "Using Variance to 
identify Important Inputs," 1994 ASA Proceedings, Section 
on Physical and Engineering Sciences, Toronto, Aug. 1994. 
S. Kaplan and L. Karklin, "Calibration of Lithography 
Simulator by Using Substitute Patterns," Proceedings on 
Optical/Laser Microlithography VI, SPIE 1927, pp. 
847-«58, 1993. 

C. Mack and E. Charrier, "Yield Modeling for Photolithog- 
raphy," Proceedings of OCG Mircolithography Seminar, pp. 
171-182, 1994. 



TMA Depict, Two-Dimensional Process Simulation Pro- 
gram for Deposition, Etching, and Photolithography, version 
3.0, Technology Modeling Associates, Inc., Palo Alto, Cali- 
fornia, 1993. \ 

"Predicting Manufacturing Variabilities for Deep micron 
Technologies: Integration of Process, Device, and Statistical 
Simulations, " in Simulation of Semiconductor Devices and 
Processes, 5, S. Selberherr, H. Stippel and E, Strasser, eds, 

Ta p. 229 -232, Springer-Veriag, New York, 1993. 
W, D. Heavlin and G.P. Finnegan, "Dual Space Algorithms 
for Designing Space-filing Experiments, " Interface 1994, 

^Research Triangle, North Carolina, Jan. 1994. 
jenichi Taguchi, "Introduction to Quality Control," p. 75 
(Asian Productivity Organization, 1986). 
Mandel, J., "The Statistical Analysis of Experimental Data," 
Chapter 12, Interscience Publishers, Wiley & Sons, New 
York, 1964. 

A. B. Owen, "Controlling Correlations in Latin Hypercube 
Samples," Journal of the American Statistical Association, 
vol. 89, No. 428, pp. 1517-1522, Dec. 1994. 

(List continued on next page.) 

Primary Examiner— ^'HM^m Grant 

Assistant Examiner — Edward F. Gain, Jr. 

(74) Attorney, Agent, or Firm — Fliesler Dubb Meyer & 

Lovejoy LLP 



(57) 



ABSTRACT 



A method of defining a toleranced process based on a 
nominal process, applicable to any manufacturing process 
wherein the output is dependent on a process having a 
number of input factors which are subject to variation, and 
have a mean and standard deviation. The method comprises 
the steps of: representing the variability of the response of a 
system to the current variability of at least one of the factors 
and at least one tightening factor; evaluating the components 
against a schedule of tightening factors in an array; model- 
ing the output of the evaluation using interpolation to 
determine a mathematical function defining the toleranced 
process; and applying at least one of the models to a 
fabrication system to complete the toleranced process. 

32 Claims, 11 Drawlag Sheets 
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STATISTICAL PROCESS WINDOW DESIGN 
METHODOLOGY 

BACKGROUNi:) OF THE INVENTION 
1. Field of Ihe Invention 

The invention is directed to the mass production of 
semiconductor devices. The invention is more particularly 
directed to the problem of perceiving interactions between 
interrelated process steps on a mass-production line and 
maintaining acceptable critical dimensions across each die 
of a plurality of mass-produced integrated circuit wafers. 

CROSS REFERENCE TO RELATED 
PUBLICATIONS 

The following publications are cited here for purposes of 
reference: 

(1) S. Kaplan and L. Karklin, "Calibration of Lithography 
Simulator by Using Substitute Patterns," Proceedings on 
Optical/Laser Mircrolitography VI, SPIE 1927, pp. 
847^58, 1993. 

(2) C. Mack and E. Charrier, "Yield Modeling for 
Photolithography/' Proceedings of OCG Microlithogra- 
phy Seminar, pp. 171-182, 1994. 

(3) TMA DEPICT, Two-Dimensional Process Simulation 
Program for Deposition, Etching, and Photolithography, 
version 3.0, Technology Modeling Associates, Inc., Palo 
Alto, Calif., 1993. 

(4) Mandel, J, The Statistical Analysis of Experimental Data, 
Chapter 12, Wiley, New York, 1964. 

(5) Z. Krivokapic and W. D. Heavlin, "Predicting Manufac- 
turing Variabilities for Deep micron Technologies: Inte- 
gration of Process, Device, and Statistical Simulations," 
in Simulation of Semiconductor Devices and Processes, 5, 
S Selberherr, H Stippel and E Strasser, eds, pp. 229-232, 
Springer- Verlag, New York, 1993. 

(6) W, D. Heavlin and G. R Finnegan, "Dual Space Algo- 
rithms for Designing Space-filling Experiments," Inter- 
face 1994, Research Triangle, North Carolina, June 1994. 

(7) B. D. Ripley, Spatial Statistics, pp. 44-75, Wiley, New 
York, 1981. 

(8) A. B. Owen, "Controlling Correlations in Latin Hyper- 
cube Samples, " Journal of the American Statistical 
Association, Vol. 89, No. 428, pp. 1517-1522, December 
1994. 

(9) W. D. Heavlin, "Variance Components and Computer 
Experiments," 199A ASA Proceedings, Section on Physi- 
cal and Engineering Sciences, Toronto, August 1994. 

(10) A. R. Neureuther and F. H. Dill, "Photoresist Modeling 
and Device Fabrication Applications," Optical and 
Acoustical Microelectronics, pp. 223-247, Polytechnic 
Press, New York, 1974. 

(11) F H. Dill, J. A. Tuttle, A. R. Neureuther, "Modeling 
Positive Photoresist," Proceedings, Kodak Microelec- 
tronics Seminar, pp. 24-31, 1974. 

(12) C. Mac, "Development of Positive Photoresists," Jour- 
nal of the Electrochemical Society, Vol. 134, January 
1987. 

(13) M. Stein, "Large Sample Properties of Simulations 
using Latin Hypercube Sampling," Technometrics, Vol. 
29, No. 2, pp. 143-151, May 1987. 

(14) M. D. McKay and R. J. Beckman, "Using Variance to 
Identify Important Inputs," 1994 ASA Proceedings, Sec- 
tion on Physical and Engineering Sciences, Toronto, 
August 1994. 

(15) William D. Heavlin and Luigi Capodieci, "Calibration 
and Computer Experiments," 1997 American Statistical 
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Association Proceedings, Section on Physical and Engi- 
neering Sciences, Anaheim, August 1997 (scheduled for 

publication Summer 1998), pp. 58-63. 

Each of the aforementioned publications is hereby incor- 
5 porated by reference. 

2. Description of the Related Art 

Modem, high-density, integrated circuit devices are typi- 
cally mass-produced with large numbers of critically- 
dimensioned features. In manufacturing, it is desirable to 
10 maintain the respective critical dimensions of each die 
within a plurality of mass-produced IC wafers constrained to 
certain respective values in order to assure desired operating 
speeds and operational characteristics of the produced IC. 

Each feature on each I C die in a mass-produced wafer is 
15 the product of a succession of many process steps. Each 
process step is controlled by a combination of variable 
process parameters. 

Different combinations of variations in process parameter 
can occur on a random basis across the numerous process 
20 steps of a mass-production line, on a die-by-die basis. This 
introduces noise into the uniformity of the product outflow 
of the production line. Sometimes a specific combination of 
process parameter deviations is relatively innocuous, some- 
times it is not. It all depends on which process parameters 
25 are deviated for a given IC die and how their respective 
process steps interrelate to establish critical dimensions on 
that given die. 

One example of such multiple, interrelated process steps 
that may be useful in understanding the interaction of factors 
30 in the present application, are those typically employed to 
define a pattern of conductive lines deposited across an 
insulator of an IC chip. 

First, a dielectric layer of a generally non-planar form is 
created across the wafer. Such a non-planar dielectric layer 
35 may constitute the combined gate oxide and field oxide of a 
CMOS device. The non-planarity of the dielectric layer may 
alternatively be attributed to the non-planarity of underlying 
trenches, steps, mesas or other topographic features of the 
chip. 

40 Next, a polysilicon or other conductive layer is deposited 
conformably on top of the dielectric layer. An anti-reflective 
coating (ARC) may be optionally deposited on the conduc- 
tive layer to reduce undesired reflections in a following 
exposure step. 

45 The deposition of the conductive layer and optional ARC 
layer is followed by a spinning-on or other deposition of a 
photoresist (PR) layer. The PR layer may or may not be 
planarized depending on process specifics. 

The photoresist-coated wafer is then positioned within a 

50 stepper by an alignment mechanism. Tiled areas of the 
photoresist (PR) are successively exposed to a stepped 
pattern of resist-modifying radiation. After the step-wise 
exposure, the photoresist is "developed" by, for example, 
soft baking the wafer to induce cross-polymerization in the 

55 photoresist material and by subsequently dissolving away 
non-hardened portions of the photoresist with a specific 
solvent. The developed photoresist defines an etch mask. 

The developed wafer is next etched, for example in a 
plasma etch chamber, so as to transfer the hardened image 

60 in the photoresist to the conductive layer. This produces a 
corresponding pattern of lines in the conductive (e.g., doped 
polysilicon) layer. 

The photoresist mask is stripped off, or kept depending on 
process specifics, and further process steps follow. One 

65 example of a further process step is the selective implant of 
dopants into exposed semiconductor regions so as to create 
self-aligned source and drain regions at opposed sides of 
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each conductive line, where the conductive line lies over 
gate oxide. The width of the conductive line at such a region 
of dopant implant defines the channel length of the formed 
IGFET transistor. 

Within each of the above-described process steps, there 
are one or more variable physical attributes (or "process 
parameters") that control the final outcome of the produced 
device. Some process parameters may be adjusted by a line 
operator. Some are "set" by the design of the process 
equipment that is installed into the mass-production line. 

Post-exposure development time and temperature are also 
subject to variance away from pre-established goal values. 
The diffusion length of the development chemistry may vary 
across a wafer. In a subsequent plasma etch, the variables 
can include: time, pressure, temperature, flow rate, and 
field-proximity effects resulting from the pitch and step 
profile of closely spaced mask features. 

Because successive steps of IC production tend to be 
interdependent, a slight variation in parameter(s) of one 
process step can be magnified by a further variation in the 
parameters of a second process step to produce unacceptable 
numbers of defective product at the output end of the 
mass-production line. 

For example, if PR (photoresist) thicknesses decreases 
slightly and the focal depth of the exposure optics also 
decreases slightly and the exposure dosage also decreases 
slightly during production of a first-sampled IC chip as 
compared to the corresponding process parameters for a 
second-sampled IC chip, the combined effect may be to 
significantly shift the position and intensity of the radiation 
exposure pattern relative to the photoresist layer during the 
production of the first-sampled IC chip. The second -samp led 
IC chip may come off the production line in acceptable form 
while the first-sampled chip comes out of the same mass- 
production line in defective form. 

It is hard to pinpoint why mass-production yield for a 
given circuit layout on a given mass-production line 
becomes unacceptably low. The statistical variance of PR 
thickness across the production lot may be small. The 
statistical variance of focus across the production lot may be 
small. But the physical interaction between the two noise 
quantities can be such that the over all lithography process 
produces chips having a much larger variance in terms of 
critical dimensions. If the overall variance in final critical 
dimensions becomes too large, production yield may suffer 
significantly. 

The interdependence of variance of a number of specific 
process steps is not easily perceived. It is not consistent 
across all steps, or even as between different product 
designs. By way of example, compare the layout of a 
densely packed memory circuit against the less regular 
layout of a random logic circuit. Some circuit layouts give 
rise to more so-called "field proximity effects" than others. 
(Field proximity effects occur in process steps such as 
exposure and plasma etch.) These field proximity effects can 
affect mass-production yield, as can many other factors. 

In short, a fairly complex matrix of relationships exists 
between material deposition operations, mechanical 
operations, optical operations, thermal operations, chemical 
reactions and the specific layout of each product. It can turn 
out for each given layout that relatively small variances in 
certain, key process parameters have far more dominant 
effects on critical dimensions of the final product than do 
larger variances in other process parameters. 

The goal of process engineers is then to determine the 
inherent trade-offs between the interdependent steps. In this 
endeavor, computer experimentation has aided the engineer 
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by allowing distributions of key output parameters to be 
estimated. Generally, this is accomplished by varying 
selected inputs or noise factors, then by running computer 
code to simulate the variation likely to occur in manufac- 
5 turing. 

Three different process analysis and design approaches 
exist with respect to the description, modeling and modifi- 
cation of a given process based upon the variation likely to 
occur in the process. The first approach, generally referred 

10 to as system design (as defined by Genichi Taguchi in 
"Introduction to Quality Control," p. 75 (Asian Productivity 
Organization, 1986)), comprises the general product 
architecture, schematic, functionahty, and feature set. In the 
second approach, generally referred to as parameter design, 

15 as defined by Taguchi, id. at p. 76, nominal process settings 
are determined in order to minimize performance variability 
and costs. In a third approach, referred to as tolerance design 
(Taguchi, id. at p. 78), tolerance and specification limits are 
determined in order to achieve engineering objectives and 

20 minimize costs. 

Computer experiments have emerged as one way engi- 
neers can model process designs using statistical methods to 
devise process parameters for semiconductor manufactur- 
ing. The models allow engineers to increase product yields 

25 by selecting variances in certain process factors such that the 
variance(s) of the output characteristic(s) does(do) not 
adversely affect product performance. 

One modeling paradigm involves a description of the 
effect of a change in one process variation on the output of 

30 another relative to some common reference. For the pur- 
poses of this description, this paradigm will be defined as an 
"engineering process window." In one aspect, a process 
window is the inherent effect of a change in one parameter 
resulting in an effect on any other parameter which may be 

35 plotted in two (or perhaps more) dimensions. Stated still 
another way, in the case of two input parameters, a process 
window is any description of the tradeoff between the 
tolerance of one input parameter and the tolerance of other 
input parameters. Parameter design gives a process window 

40 which defines, as a function of input parameters Xj . . . Xy> 
the region in which the product performs in conformance to 
specifications. Tolerance design yields, for a given range of 
input parameters, the tolerances X^±A^, . . . , X^Ay which 
allow the product to perform in conformance to its perfor- 

45 mance specification. 

Various process windows may be defined by perturbing 
input factors by a fixed amount and determining the average 
desirability of response values over a grid of points, such as 
the number of variables which exist in a semiconductor 

50 fabrication process. 

FIG. 1 illustrates an exemplary engineering process win- 
dow. In FIG. 1, three linewidths 0.2925 fiva, 0.325 /on, and 
0.3575 fim resulting from a particular process are plotted 
based upon the required values of photoresist exposure 

55 (y-axis) and focus (x-axis) necessary to achieve each lin- 
ewidth. Here, exposure times are defined in seconds, and 
focus in micrometers. In this figure, the "goal" linewidth is 
0.325 /im, while the acceptable variation is 10% (or 
±0,0325) yielding the narrow (0.2925) and wide (0.3575) 

60 linewidth limits. In this instance, a process window 30 is 
defined as an acceptable range of inputs of focus and 
exposure (as defined relative to the x-y axes) which yield a 
linewidth within this 10% range. The tolerance is the varia- 
tion of one factor which yields the acceptable output within 

65 the window 30 relative to the other component's value. 
In computer experimentation involving engineering pro- 
cess windows, calibration to realistic measurements is gen- 
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erally compared with computer implemented models. Data 
from the process may be sampled from sacrificial areas of 
actual processed substrates, with the sample data then used 
with empirical experiments or computer simulators (such as, 
for example, the DEPICT photoHthography simulator avail- 
able from Technology Modeling Associates of Palo Alto, 
Calif.) to provide a calibrated model of the process results. 

FIG. 2 represents a contour plot comparing three different 
foci from a calibrated computer experiment (CAL), an 
un-calibrated exposure/focus matrix of a simulation 
(UNCAL), and measured experimental results, of the expo- 
sure latitude versus the depth of focus in a line -width 
definition process. As will be noted therein, computer 
experimentation both calibrated process simulations and 
computer experimental methods have a relatively close 
accuracy to actual results. 

Process windows which compare results of altering two 
factors are relatively simple; the model deals with changes 
to one parameter which affect another parameter. The analy- 
sis becomes much more complex where a number of input 
factors are involved as in, for example, a cross-correlation of 
device parameters which are utilized in a semiconductor 
fabrication process and affect conductive line formation. 
Such factors include, and are abbreviated in this specifica- 
tion as follows: 
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where X=0, perfect control (no variance) of the input factor 
is presumed. One factor per iteration is presumed perfect, 
and all other factors are allowed to vary. 

FIG. 3B is a representation of the resultmg quantitative 
AMV analysis for seven factors: exposure, focus, MaskCD 
(the exposure mask thickness critical dimensions), prThk, 
BARC, PEB, and devT for an exemplary sample. In FIG. 
3B, Syis the LHC-based estimate of the variation attributable 
to factor f, Sq^ the estimate of the total variation induced 



10 from all factors, and S 



15 



20 



25 



35 



45 



One descriptive method used successfully to compute 
variance parameters is known as analysis of manufacturing 
variance (AMV) and is set forth in Heavlin, "Variance 50 
Components and Computer Experiments," reference paper 
no. 9, cited above. 

In analysis of manufacturing variance, the goal is to 
decompose the total variance of a process into components, 
one component associated with each input factor, and con- 55 
sistent with the variation to be anticipated in manufacturing. 
In a basic form of AMV, each noise factor is perfectly 
controlled, and the reduction in output variation noted. AMV 
allows for control of subsets of factors as well. 

AMV in its simplest form comprises using Latin hyper- 60 
cube (LHC) samphng for factors . . . in a computer 
simulation, calculating the total variance of the resulting 
values, and substituting the mean of each factor to calculate 
the reduction in variance from the total variance to estimate 
the contributions of the factors held constant. AMV analysis 65 
is illustrated in the table shown in FIG. 3A. In FIG. 3 A, each 
>^ represents a noise factor for an input under determination 



^the component of the total 



variance that can be attributed to factor f. As shown in FIG. 
3B, the AMV analysis results in a showing that the three 
most important factors contributing to the total variation in 
the process under consideration are the post exposure bake 
time (PEB), exposure time (Exp) and focus. Given the 
results of FIG. 3B, the conclusion drawn for the process 
considered therein is that the most attention should be given 
to controlling the post exposure bake time and exposure 
dosage variance in the corresponding manufacturing process 
under consideration. By way of example, driving exposure 
dosage variance down to zero provides a 38% improvement 
for the overall cross-reticle variance of a dense-line circuit 
layout. 

AMV does, however, exaggerate the contributions from 
third order and higher order interaction terms. FIG. 4 
illustrates a calculation of the bias of AMV's interaction 
term (i.e. the tenn being held constant). In FIG. 4, dots 
represent "l'"s and blanks, "0"'s; o^-o^ represent the 
standard deviation of the interaction terms; and T^-T^23A5 
the iteration for each factor. 

AMV nevertheless provides a strong methodology for 
deriving a quantitative description of the interaction between 
a number of factors. AMV does not provide an effective 
method for deriving the tolerances of each input factor 
which are suitable for maintaining the operating character- 
istics and yield of the final product. AMV quantifies the 
leading sources of variation, relative to the control assump- 
tions on the noise factors. Based on computer experiments, 
applying AMV requires one to consider the costs of running 
the computer code, the non-linearities of the variance 
structure, and the number of factors under evaluation. 

Another methodology related to AMV is described in U.S. 
Pat. No 5,646,870 and referred to therein as incremental 
leveraging. Initially, a predefined goal for variation is 
selected — for example 10% — and an evaluation of each 
factor made on this basis. Process parameters are selected 
one at a time and the variance of the selected process 
parameter is temporarily reset to zero. The noisy set of 
simulations is re-run and the new total variance is calculated 
and stored while still retaining the variance value of an 
original first run where none of the relevant variances were 
reset to zero. 

FIGS. 5A and 5B illustrate the results of such an incre- 
mental leveraging process. As seen in FIG. 5A, the goal of 
the operation is to reduce the 3 a/mean value for the overall 
process from 0.20 to 0.10 (or 10% (logarithmic). In both 
FIGS. 5 A and 5B, the x-axis increments represent a 10% 
incremental tightening of each input. In FIG. 5 A, the 
3a/mean value on the y-axis represents the response of the 
process under consideration to the tightening of the factors. 

In FIG. 5B, the first iteration of the leveraging algorithm 
pinpoints the variance of the focus as the tightening factor 
which provides the most beneficial improvement. Develop- 
ment time and mask CD remain constant over the first 5 
increments (measured along horizontal axis). At the sixth 
increment, development time provides the most beneficial 
improvement. Best beneficial improvement from mask criti- 
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cal dimension variance does not arise until the 15lh incre- function to determine suitable tolerance models impacting 

ment. The goal for the total variance of the production line control of one or more of said input factors. The method also 

is reached at the 12th increment, however, by flip-flopping includes defining a toleranced process by applying at least 

tightenings the development time and the focus parameters. one of the models to the manufacturing process steps by 

At every increment along the x-axis, each noise factor is 5 setting the value of said components and manufacturing the 

lightened in turn by ten percent. That noise factor giving the device using the toleranced process, 

greatestreductionof CD variation is selected for permanent More generally, the method of the invention can be 

tightening. The predicted critical dimension from such tight- applied to any manufacturing process wherein the output is 

ening is plotted on the upper panel of HG. 5A; the lightening dependent on a process having a number of input factors 

multipher for that noise factor is plotted on the lower panel lO which are subject to variation. In particular, a method of 

of FIG. 5B. defining a toleranced process based on a nominal process 

In the illustration HGS. 5A and 5B, incremental lever- having a plurality of components, each component having a 

aging is performed over 20 improvement increments. For mean and standard deviation; comprising the steps of: 

the first five increments, the factor with greatest leverage is representing the variability of the response of a system to the 

focus, for the next two, development time, and so on. In is current variability of at least one of said factors and at least 

total, focus is assigned 7 improvement increments and one tightening factor; evaluating said components against a 

development time 5 before we achieve 10 percent (3 a) schedule of tightening factors in an array; modeling the 

linewidth control. We achieve this level by lightening focus output of the evaluation using interpolation to determine a 

by a factor of 0.48 (0.9'') and development time by a factor mathematical function defining the toleranced process; and 

of 0.59 (0.9 ). 20 applying at least one of the models to a fabrication system 

Although the incremental leveraging techniques of FIGS. to complete the model. 
5A and 5B are applied to variance lightening, a similar 

approach may be taken to selecting the best goal settings for BRIEF DESCRIPTION OF THE DRAWINGS 
the mean process parameters such as focus depth and 

photoresist thickness. Thus, if something goes wrong on the 25 invention wiU be described with respect to the 

production line after new equipment has been installed, line particular embodiments thereof. Other objects, features, and 

operators can use the results of the incremental leveraging advantages of the invention will become apparent with 

algorithm to determine which process parameter requires reference to the specification and drawings m which: 

tweaking first and by how much. FIG. 1 is a plot representing the interaction between the 

Incremental leveraging has the advantage that it is good at 30 factors of exposure and focus in a process defining a 

suggesting one solution for input tolerance design. It does semiconductor device linewidth relative to three actual 

not, however, provide a mechanism for flexibility in allow- linewidths, 

ing more than the one solution it provides for changing the FIG. 2 is a contour plot of the exposure latitude versus 

tolerance constraints of the factors in the process under depth of focus comparing a computer experimental model, 

consideration. 35 actual measured results, and a computer simulation. 

cTT\rf\^ADvncTnjc TKn^ivmrMvT ^ ^° iUustrative table representing the factor 

SUMMARY OF THE INVENTION ^^^^y^^^ ^^^^^^ ^^^^^^^ manufacturing variance meth- 

Hence, an objective of the invention is to allow process odology. 

engineers to determine tolerances for a multitude of factors. FIG. 3B is a table representing the output of an analysis 

A further objective of the present invention is to provide of manufacturing variance analysis on a given process, 

an evaluation method for process engineers to make quali- FIG. 4 is a table illustrating the effect of higher order 

tative judgements about the manufactur ability of certain interaction terms in an analysis of manufacturing variance 

processes based on models of the output given a variety of analysis. 

inputs to the process. ^5 FIGS. 5A and SB illustrate incremental leveraging algo- 

Yet another object of the invention is to improve yields in rithm analysis of a manufacturing process. 

multi-faceted manufacturing processes, and particularly piGS. 6A and 6B illustrate the input, and output, 

semiconductor manufacturing applications. respectively, of an incremental averaging algorithm process 

In one aspect, the invention, roughly described, comprises over 210 intervals. 

a method of fabricating a mulli^omponent semiconduct or 50 pic. 7 is a flow chart depicting one embodiment of the 

device assembly. This method is composed of the steps of method of the present invention. 

first defining a nominal semiconductor manufacturing pro- ^^1- c • • * ^* 1 * ui n * *• 

™ . ^ . , . 1 I . r f FIG. 8 is an mteractive dialog table illustrating a portion 

cess. This nominal process has a plurahly of process steps r,, *u j r *u * • . r *if • i-t^^ o 

J J • i.1 ■ 1 of the method of the present invention as set forth m FIG. 8. 
used to produce the device assembly. The nominal process 

also includes a number of process input factors X^-X. 55 "^^^^ ^ ^ ^^art illustrating the statistical process 

which, when perfectly controlled to meet the nominal pro- wmdow evaluation portion of the method of the present 

cess sellings for such factors, produce a device assembly uivention. 

with defined operating characteristics. The input factors FIG- 10 is a graphical representation of the command 

each have a distribution of acceptable tolerances which are instructions and matrix input for the method set forth in FIG. 

acceptable inputs to the process and which result in the 60 9. 
defined operating characteristics. The method further 
includes evaluating the nominal process by representing the 
variability of the response of the operating characteristic as 

a function of the current variability of at least one of said In manufacturing processes, such as those utilized for 

factors and at least one lightening factor, evaluating said 65 manufacturing a semiconductor device, imperfect control of 

input factors subject to a variety of lightening factors array input variables propagates through the manufacturing pro- 

over [0,1/, and modehng the output using an interpolation cess to induce variation in the manufactured device, and in 
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response variables describing the device. Where there are a 
number of process inputs, one can define the relative impor- 
tance of such inputs, and more than one set of input 
tolerances which will meet the desired response. This means 
that there are trade-ofik in the control tolerances, or required 
distributions of standard deviations, among the various input 
factors. While process windows represent one-to-one trade- 
offs, a more generalized model is required to consider all 
sets of input tolerances. 

While the method of the present invention will be 
described herein with respect to a semiconductor manufac- 
turing process, it should be recognized that the method of the 
present invention may be applied to any process wherein a 
number of process variables affect the output variance of the 
manufactured product. 

In accordance with the present invention, incremental 
leveraging methodology is utilized as a basis to define a 
statistically based, process window algorithm wherein the 
response variability of a process is represented as a function 
of the tolerances of each variable. In brief, one utilizes the 
tolerance variations defined in the aforementioned algo- 
rithms as inputs to determine, for any number of process 
variables, the effect of varying the tolerance of each factor 
on the tolerances of other factors under consideration rela- 
tive to a response factor for the entire process. This algo- 
rithm may be performed in a computer and have an output 
provided to a process engineer who may then determine 
optimal modifications to the nominal process under consid- 
eration to derive a toleranced process. 

In general, the method of the present invention comprises: 
defining a nominal process; evaluating F input factors 
{X^f=l , . . . , F} which affect the output of the nominal 
process; and defining, as a result of the evaluation, a 
toleranced process having an output response less suscep- 
tible to variance in the input factors. In one aspect, the 
invention comprises a method for manufacturing a semicon- 
ductor device which comprises: defining a nominal semi- 
conductor manufacturing process; evaluating a subset of 
factors controlling process steps in the manufacturing pro- 
cess; modifying the process to reduce the effect of large - 
effect factors by constraining the control on such factors 
and/or loosening control small-effect factors; and manufac- 
turing a semiconductor device based on the modified pro- 
cess. 

In brief, the evaluation portion of this method comprises: 
representing the variability of the response of the process as 
a function of the factors' current variability and selected 
tightening factors; performing an evaluation where the tight- 
ening factors are varied using an orthogonal array, varying 
the factors over [0,1]''; and modeling, using kriging 
interpolation, a neural network or equivalents, to determine 
moments or other statistics of a distribution of the process 
for given changes in the constraints of each factor. 

The method of forming the semiconductor device in 
accordance with the present invention results in better 
planning of next-generation equipment requirements and 
ultimately in improved product yields as the considered 
factors are optimized for the given process by allowing 
modification of the nominal process factor inputs based on 
the model resulting from the evaluation step. 

By way of example, a typical semiconductor mass- 
production line which is part of an overall wafer processing 
system may include: (1) a station for coating each wafer with 
material layers such as the anti-reflective coating and pho- 
toresist layers; (2) a station for exposing the material layers 
(the photoresist layer) to a radiation pattern; (3) a station for 
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developing each exposed wafer; and (4) a station for etching 
each developed wafer. The method of the present invention 
allows tweaking of the process at any one or multiples of 
such stations. 

5 Typically, each post-etch wafer lot is moved further 
downstream for further processing (e.g., dopant implant) in 
accordance with well known techniques. Before such further 
processing, sample wafers from each post-etch wafer lot 
may be sampled to determine critical dimensions in order to 
calibrate the computer model of the method of the present 
invention. 

The evaluation method of the present invention may be 
performed in an independent computer unit or can be a 
module within a computer that also carries out functions of 

J 5 other described or to-be described modules. The output of 
the model may be used by the process engineer in possession 
of other data, such as cost factors, to define constraints and 
make choices for altering the semiconductor manufacturing 
process to determine which constraints to confine or relax. 

20 According to the evaluation portion of the method of the 
invention, a statistical process window evaluation is used 
over a set of F individual factors. For each component, the 
simulation is used to generate a statistical distribution of 
each fit-determining parameter. The distributions have a 

25 mean and a deviation about the mean, and also between 
manufacluring-capability-defined lower and upper limits, 
e.g., mean ±3 standard deviations. 

FIGS. 6A and 6B illustrate the output of the incremental 
leveraging calculation for a given simulation. The input 

30 table 40 includes the nominal mean and the standard devia- 
tion values for seven input factors considered in this 
example which comprises line width definition process, 
similar to that described above and in U.S. Pat. No. 5,646, 
870. Each factor is assigned an effort of 0.9, reflecting that, 

35 for this exercise, the design goal is a 10% incremental 
change. It should be recognized that other effort factors may 
be utilized in accordance with the present invention. The 
factor of a 10% increment is selected based upon common 
practice in the industry. The output table 50 shown in FIG. 

40 6B comprises a calculation in accordance with the incre- 
mental leveraging algorithm described with respect to FIGS. 
5A-5B, over 210 increments. Column 1 of output table 50 
shows the resulting weighting of weighted tolerance factor 
for each of the seven factors analyzed. For example, for 

45 incremental row 205, the focus factor has a standard devia- 
tion of 0.06960, and the factors of exposure (ExpD), focus 
and post exposure bake (PEB) are shown to have the greatest 
contribution over these increments. Incremental leveraging 
assigns the focus parameter 7 increments, exposure 11 

50 increments, and PEB 12 increments. This yields tightening 
factors for ExpD of 0.3138 1-(0.9)^\ for Focus of 0.47830= 
(0.9)^ and for PEB of 0.28243 =(0.9)' ^ This output of the 
incremental leveraging algorithm provides the tightening 
factor starting points for determining the matrix of tighten - 

55 ing factors used in the evaluation portion of the present 
invention. 

FIG. 7 is a flowchart depicting, in one aspect, a method in 
accordance with the present invention. At step 55, a manu- 
facturing process under consideration is defined based on the 

60 particular characteristics of the product to be manufactured. 
In the case of a semiconductor device, this process will 
include a number of process steps as described herein. The 
process defined to produce the product is generaUy referred 
to as the "nominal" process and results in process steps, each 

65 process step having one or more input factors, such as those 
factors set forth above with respect to the semiconductor 
manufacturing process. 
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As shown at step 60, for a set of process steps, the input 
factors whose variance affects the output yield of the process 
based on their interaction with other factors are defined. 
Each input factor will have a nominal or mean value, and a 
standard deviation about the mean value. The standard 
deviation can be an actual deviation during the process or a 
simulated deviation for purposes of the evaluation of the 
present invention. 

Step 100, defined in dashed lines and comprising several 
sub-steps, comprises the evaluation portion of the method of 
the present invention. The evaluation portion 100 can com- 
prise a set of instructions which are processed by a 
computer, in one or more modules as described above, 
designed specifically to implement the evaluation based on 
the input factors provided to it. 

In one embodiment, as depicted in FIG. 7, the evaluation 
portion 100 of the method of the present invention 
comprises, at step 110, defining the variance of the manu- 
facturing output — the product produced — as a function of 
the variability of the input factors and tightening factors 
applied to each input factor. Once so defined, at step 120, an 
evaluation is performed by varying the tightening factors to 
determine alternative input tolerances which yield accept- 
able manufacturing outputs. Finally, at step 130, the evalu- 
ation is modeled to determine the process variance as a 
function of the tightening factors. 

Once the evaluation is performed, at step 70, the nominal 
manufacturing process may be varied lo derive a toleranced 
manufacturing process by implementing a tolerance model. 
Finally, the toleranced process is used to fabricate a device 
at step 80. 

FIG. 8 shows an interactive dialog for a statistical process 
window algorithm that is used in the evaluation portion of 
the present invention. FIG. 9 presents a generalized descrip- 
tion of the evaluation algorithm of the present invention 
which should aid in understanding the more specific descrip- 
tion to follow. In general, the dialog depicts the interaction 
between terms necessary to loosen the tightening factor for 
focus from an interactive leveraging-derived value to 1.0, 
and the effect on the tightening factors of all other factors 
under evaluation. It should be noted that while focus is 
illustrated as the factor under consideration, other factors, or 
multiple factors, may be considered. 

The tightening factors for exposure, focus and post 
expose bake are calculated tightening factors based on an 
incremental leveraging algorithm, as described above. All 
other factors are assigned a factor of 1.0=(0.9)°. A resuhing 
control factor is computed. The control factor reflects the 
3a/mean value as computed in an incremental averaging 
algorithm. Next, the tightening factor for focus is relaxed 
(from 0.47830 to 0.081 (at coL 2)). and all other factors 
except ExpD and FEB tightened, and the resulting control 
recalculated. As will be noted, the control factor increases to 
0.11536. Next all factors, including ExpD and FEB are 
tightened, and the control factor again calculated. This 
process continues as illustrated in FIG. 8. 

The goal of the aforementioned dialog can be seen by 
reviewing the subsequent rows 4-11 wherein the focus 
contribution is relaxed to zero increments (0.9)^«1 .0 in row 
6, as all other factors are tightened until the resulting control 
factor achieves proximity with the initial, incrementally- 
derived control factor in row 1 of the dialog. 

Comparing the statistical process window dialog with 
incremental leveraging, incremental leveraging gives an 
output for each factor of (o) (lY for each tolerance. With 
SPW, the resulting control allows one to relax or modify any 
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number of the factors (not merely the single focus factor 
identified in the dialog of FIG. 8), to achieve the resulting 
control value. It should be noted that in the method of the 
present invention, each of the process factors may be varied 
5 relative to external considerations, such as cost or ease of 
control over the factor. In the aforementioned example of 
FIG. 8, the noise factor is selected at (0.9). In the method of 
the present invention, the noise factor is varied over [0,1]^ 
(mathematically, the interval from 0 to 1, inclusive) in an 
orthogonal array by means of a computer simulation as 
described below. 

The evaluation portion 100 of the method of the present 
invention will be described with reference to FIGS. 9 and 10. 
FIG. 9 is a flow chart of the evaluation algorithm 100 of the 
present invention FIG. 10 is a representation of the com- 
mand instmctions and mathematical inputs of the evaluation 
method of the present invention. 

FIGS. 9 and 10 consider, for example, a lithography 
portion of a semiconductor manufacturing process. A no mi- 
nal lithography process has been determined by process 
engineers, based on the device to be manufactured, and the 
nominal process will be evaluated in accord with the inven- 
tive method. 

As shown in FIG. 9, process step UOa consists of defining 

25 a resulting linewidth critical dimension as a mathematical 
function of a number (f) of input factors. The result is an 
output variable (cd„dc5a) which is the goal of the linewidth 
definition process (e.g., a linewidth). In one aspect of the 
present invention, this definition step 11 Oa can be performed 

3Q by a kriging interpolation in accordance with the approach 
outlined in B. D. Ripley, "Spatial Statistics," article 7 set 
forth above, or by suitable fast computer simulation tech- 
niques. Essentially, step llOa can be decomposed into the 
substeps of: building a kriging model to determine the 

35 response dimension relative to the input factors; and cali- 
brating the model. The calibration adjustment may be per- 
formed in accordance with the teachings of W. D. Heavlin 
and L. Capodieci in "Calibration and Computer 
Experiments," cited as article 15, above, wherein simulated 

4Q critical dimensions are compared against empirical critical 
dimensions. The calibration model may result from either 
kriging inteipolation or use of a neural network in accor- 
dance with the teachings of the aforementioned reference. 
Hie output critical dimension used in the method of the 

45 present invention (cd_dc5fl) is the result of the calibration 
model and the simulation critical dimension resulting from 
the kriging model. 

In FIG. 10, the evaluation steps llOa^ 115 and 120a are 
represented as a single command fine 90 which may be used 

50 to call the functions or tables described herein. The tables 
and arrays described with respect to FIG. 9 are graphically 
represented. In FIG. 9, at step 115, the nominal process 
values of the input factors and the starting process tolerances 
are defined for input to the evaluation portion of the method. 

55 The input factors (mean) and variances (stdev) are set forth 
in a table (pw_mst). The nominal process targets at the input 
factors are the staning (wide-standard deviation, current 
generation) process tolerances of the input factors. 
The current generation values of the input factors are then 

60 used to build a data set for the evaluation portion of the 
method of the present invention. 

At step 120^^ a schedule of tightening factors (L49 as 
illustrated in FIG. 10) is then evaluated against the defined 
output variable (cd_dc5fl), which is subjected to a Latin 

65 hypercube sampling (LHS343) for each row of the schedule 
of tightening factors in accordance with the following 
pseudocode: 
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loop "schedule" (L49) What is claimed is: 

for row i of schedule 1- A method of fabricating a multi-component semicon- 

copy LHS,. from LHS ductor device assembly, comprising: 

loop j (A) defining a nominal semiconductor manufacturing 

Col j of LHS-mst(mean)Q)+(Col j of LHS,-)*mst 5 process having a plurality of process steps which 

(stdXj)*sched(i,j) produce the device assembly, the nominal process 

end loop j * including a number of process input factors X^-X^- 

I H<; nA-i. ro«/c^^r-H Hf-^/i which, when perfectly controlled to meet the nominal 

i.rtj>,p^jrowsji ca_ac3fl ^^^^ settings for such factors, produce a device 

m(i)=mean of UiS343-od dcSa assembly with defmed operating characteristics, the 

s(i)«stdev of LHS343-*cd_dc5fl 10 input factors having each having a distribution of 

end loop 1 acceptable tolerances which are acceptable inputs to the 

Step 128 results in an output (pw__dc5ci/ as represented in process which result in the defined operating charac- 

FIG. 10), which is a mean and standard deviation table with teristics* 

a number of rows (49) coiresponding to the number of rows (g) evaluating the nominal process by: 

(49) in the schedule (L49) of tightenmg factors. The output representing the variability of the response of the 

pw dcScrf IS then used m a kriging mterpolation at step characteristic as a function of the cunent 

130a to yield a mathematical function which is the standard variabiliTy of at least one of said factors and at least 

deviation of the output response and the tightening factors. . , / . - 

The output of the kriging interpolation shall be referred to igtitenmg lactor; 

herein as the "statistical process windows (SPW) function." (2) evaluating said input factors subject to a variety of 

The SPW function is the standard deviation of the output 20 tightening factore where the tightenmg factors are 

response as a function of each input factor and a correspond- varied over [0,1/ using an array; 

ing tightening factor for the particular input factor P) modeUng the output using an interpolation function 

As represented in FIG. 10, at step 94, the log of col. 1 to determine suitable tolerance models impacting 

(stdev) of py/__6c5cd is taken to exaggerate small values of control of one or more of said input factors; 

the output variance. As shown in FIG. 10, in the particular (C) defining a toleranced process by applying at least one 

kriging estimation which is utilized in the present 25 Qf the models to the manufacturing process steps by 

invention — kpw_dc5ci/, two iterations are used with a ridge setting the value of said components; and 

(a.k.a. nugget) parameter of 0.01. (D) manufacturing the device using the toleranced pro- 

At step 98 of FIG. 10, the exponent taken in accordance cess, 

with conventional statistical techniques to invert the log 2. '^The method of claim 1 wherein said step (B)(1) 

function of step 94. 3q comprises the sub-steps of: 

As shown at step 98, the kriging interpolated estimate is defining an output response based on a matrix of said 

a deviation of 0.060705, compared to a directly calculated factors* and 

value of 0.060533 for one particul^ point of interest^ calibrating said output response based on empirical data 

An additional mode of method of the invention is shown ^he outout response 

in FIG. 9 wherein, at Une 135a a process engineer may enler 3 ^^.^^j „f ^^^-^^ j wherein said sub-step of 

into an interactive mode with the method of the present 35 ^^^^^ comprises utilizing a kriging interpolation, 

invention to propose alternative schedules of tightening 4 j^,^ method of claim 2 wheTein said sub-step of 

factors based on external factors such as the cost of con- defining comprises utilizing a neural network, 

trolling a particu ar factor or other considerations and 5 ^^^^ ^g,,^^,^ ^ ^^.^^j^ ^-^ ^ (b)(2) 

determine the resultmg output based on this proposed sched- . r \ /v / 

ule of tightening factors. compr^s. ^ . ^ . ^ , . . 

In stiU another aspect of the invention, AMV or the for the array of tightening factors, evaluating for each row 

incremental leveraging algorithm can be used with the of the array, a Latin hypercube sample of, for each 

statistical process windows function. The result is a faster column of the sample, the mean and standard deviation 

evaluation than using AMV alone, but one which is less of each input factor, taking the mean of the Latin 

accurate than the above proposed methodology. hypercube sample output of each row of the array in the 

The many feaUires and advantages of the present inven- 45 standard deviation of the Latin hypercube sample for 

tion will be readily apparent to one of average skill in the art. ^^^h row of the array. 

The method of the present invention provides a model which 6- The method of claim 1 wherein said step (B)(3) 

can be used in processing devices based upon computer comprises applying a kriging interpolation to the output of 

simulated results of processes, or on empirically derived said step (BX2). 

models of each process. No specification of a cost-of- 50 ^he method of claim 1 wherein said step (B)(3) 

tolerance function is required, thereby simplifying applica- comprises utilizing a neural network on the output of step 
tion of the method to various processes or sub-processes in 

(B)(2). 

device manufacturing. 8* The method of claim 1 wherein steps (B)(2) and (B)(3) 

The method of the present invention utilizes an evaluation are repeated based on empirical choices made according to 

process which makes explicit the trade-off in the control external factors of the nominal process, 

requirements amongst various manufacturing steps, and ^ The method of claim 1 wherein at least one of the 

provides a quantitative analysis of this trade-off for use by tightening factors is selected by a user, 

process engineers in developing the fabrication process. 10. The method of claim 9 wherein steps (B)(2) and (B)(3) 

Process engineers are provided with an explicit role for are repeated, wherein at least one of the tightening factors is 

engineering judgment by the quantitative description of each selected by the user. 

trade-off in the process. ^° H- The method of claim 1 wherein said step (B)(1) 

Obvious modifications to the method of the invention will comprises: 

be apparent to those of average skill in the art. For example simulating an output response as a function of the input 

and without limitation, at step 130a, other methods of factors to determine an output variable; 

modeling the output of step I20a, such as utilizing a neural calibrating the simulation by comparing empirically 

network, are contemplated. All such modifications are 65 determined output responses as a function of selected 

intended to be within the scope of the invention as defined factors to simulated output responses generated by the 

by the claims. simulation. 
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12. The method of claim 11 wherein the step of simulating 
comprises using a kriging interpolation. 

13. A method of fabricating a multi-component assembly, 
the method including at least components Xl-Xf wherein 
each of said components Xl-Xf has a distribution of accept- 5 
able tolerances within which the method will produce the 
multi-component assembly in an acceptable fashion, com- 
prising: 

(A) defining the nominal values of each component; 

(B) evaluating the method by: -jq 

(1) defining the variance of at least one critical dimen- 
sion of the assembly as a function of the variance of 
the components and a plurality of tightening factors, 
one tightening factor associated with each compo- 
nent; 

(2) analyzing the nominal components by comparing a 
Latin hypercube sampling of the nominal compo- 
nents against a schedule of tightening factors; 

(3) modeling the analyzed components to define the 
variance of the output as a function of the compo- 
nents and the tightening factors; and 

(C) manufacturing the multi-component assembly by per- 
forming the method and controlling ones of said com- 
ponents based on the output of step (B). 

14. The method of claim 13 wherein said step (B)(1) 
comprises the sub-steps of: 

defining an output response based on a matrix of said 

components; and 
calibrating said output response based on empirical data 

of the output respon^. 

15. The method of claim 14 wherein said sub-step of 
defining comprises utihzing a kriging interpolation. 

16. The method of claim 14 wherein said sub-step of 
defining comprises utihzing a neural network. 

17. The method of claim 13 wherein said step (B)(2) 
comprises: 

for the schedule of tightening factors, evaluating for each 
row of the schedule, a Latin hypercube sample of, for 
each column of the sample, the mean and standard 
deviation of each component, taking the mean of the 
Latin hypercube sample output of each row of the 40 
schedule in the standard deviation of the Latin hyper- 
cube sample for each row of the schedule. 

18. The method of claim 13 wherein said step (B)(3) 
comprises applying a kriging interpolation to the output of 
said step (B)(2). 45 

19. The method of claim 13 wherein said step (B)(3) 
comprises utilizing a neural network on the output of step 
(B)(2). 

20. The method of claim 13 wherein steps (B)(2) and 
(B)(3) are repeated based on empirical choices made accord- 
ing to external factors. 

2L The method of claim 13 wherein said step (B)(1) 
comprises: 

simulating an output response as a function of the com- 
ponents to determine an output variable; 

calibrating the simulation by comparing empirically 
determined output responses as a function of selected 
components to simulated output responses generated 
by the simulation. 

22. llie method of claim 21 wherein the step of simulating 
comprises using a kriging interpolation, 

23. A method of defining a toleranced process based on a 
nominal process having a plurality of components, each 
component having a mean and standard deviation, compris- 
ing the steps of: 

(A) representing the variability of the response of a 65 
system to the current variability of at least one of said 
components and at least one tightening factor; 
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(B) evaluating said components against a schedule of 
tightening factors in an orthogonal array; 

(C) modeling the output of the evaluation using interpo- 
lation to determine models defining the toleranced 
process; and 

(D) applying at least one of the models to a fabrication 
system. 

24. The method of claim 23 wherein said step (A) 
comprises: 

defining an output response based on a matrix of said 

components; and 
calibrating said output response based on empirical data 

of the output response. 

25. The method of claim 24 wherein said step of defining 
comprises utilizing a kriging interpolation. 

26. The method of claim 24 wherein said step of defining 
comprises utilizing a neural network. 

27. The method of claim 23 wherein said step (B) 
comprises: 

for the schedule of tightening factors, evaluating for each 
row of the schedule, a Latin hypercube sample of, for 
each column of the sample, the mean and standard 
deviation of each component, taking the mean of the 
Latin hypercube sample output of each row of the 
schedule in the standard deviation of the Latin hyper- 
cube sample for each row of the schedule. 

28. The method of claim 23 wherein said step (C) 
comprises applying a kriging interpolation to the output of 
said step (B), 

29. The method of claim 23 wherein said step (C) 
comprises utilizing a neural network on the output of step 
(B). 

30. The method of claim 23 wherein steps (B) and (C) are 
repeated based on empirical choices made by external 
factors. 

31. The method of claim 23 wherein said step (A) 
comprises: 

simulating an output response as a function of the com- 
ponents to determine an output variable; 

calibrating the simulation by comparing empirically 
determined output responses as a function of selected 
components to simulated output responses generated 
by the simulation. 

32. A method, comprising: 

(A) defining a product and a nominal process comprising 
a series of process steps, the nominal process including 
input factors for the steps suitable for producing the 
product; 

(B) defining a set of the input factors whose variance 
affects the output characteristics of the process for a set 
of process steps; 

(C) defining the variance of the manufacturing output of 
the nominal process as a function of the variability of 
the input factors and tightening factors; 

(D) evaluating the input factors by applying a schedule of 
tightening factors to determine the tolerances of the 
input factors relative to the output factor; 

(E) modeling the evaluation to determine process variance 
as a function of the tightening factors; 

(F) defining a toleranced process for producing the prod- 
uct; and 

(G) producing the product by performing the toleranced 
process. 
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